library(dplyr) ; library(tidyr) ; library(ggplot2)root <- rprojroot::has_file(".git/index")datadir = root$find_file("data")funsdir = root$find_file("functions")savingdir = root$find_file("saved_files")files_vec <-list.files(funsdir)for( i in1:length(files_vec)){source(root$find_file(paste0(funsdir,'/',files_vec[i])))}# dat_tax = data.table::fread('https://raw.githubusercontent.com/rafaelcatoia/zoop_16N/main/treated_taxonomy_dat.csv') %>%# as_tibble()# # ### Loading the new data # grump_version = 'new'# # if(grump_version=='new'){# ## Use this if you are using the new GRUMP Data Set# datapath = root$find_file(paste0(datadir,'/','grump_asv_long20240501.csv'))# dframe = data.table::fread(input = datapath) %>%# filter(Cruise %in% c('P16N','P16S')) %>% # #filter(Raw.Sequence.Counts>0) %>% # filter(Domain!='Unassigned') %>% # mutate(Raw.Sequence.Counts = Corrected_sequence_counts)# # }else{# ## Or this one if you are using the OLD GRUMP Data Set# datapath = root$find_file(paste0(datadir,'/','grump_asv_long_20240110.csv'))# dframe = data.table::fread(input = datapath) %>%# filter(Cruise %in% c('P16N','P16S')) %>% # #filter(Raw.Sequence.Counts>0) %>% # filter(Domain!='Unassigned')# }# # dat_tax %>% select(starts_with('P16')) %>% dim()